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Abstract 

A N-sized inertial classical Heisenberg ferromagnet, which consists in a 
modification of the well-known standard model, where the spins are re- 
placed by classical rotators, is studied in the limit of infinite-range interac- 
tions. The usual canonical-ensemble mean-field solution of the inertial clas- 
sical n-vector ferromagnet (for which n = 3 recovers the particular Heisen- 
berg model considered herein) is briefly reviewed, showing the well-known 
second-order phase transition. This Heisenberg model is studied numer- 
ically within the microcanonical ensemble, through molecular dynamics. 
In what concerns the caloric curve, it is shown that, far from criticality, 
the kinetic temperature obtained at the long-time-limit microcanonical- 
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ensemble simulation recovers well the equilibrium canonical-ensemble esti- 
mate, whereas, close to criticality, a discrepancy (presumably due to finite- 
size effects) is found. The time evolution of the kinetic temperature indi- 
cates that a basin of attraction exists for the initial conditions for which the 
system evolves into a metastable state, whose duration diverges as N — > oo, 
before attaining the terminal thermal equilibrium. Such a metastable state 
is observed for a whole range of energies, which starts right below criticality 
and extends up to very high energies (in fact, the gap between the kinetic 
temperatures associated with the metastable and the terminal equilibrium 
states is expected to disappear only as one approaches infinite energy). To 
our knowledge, this has never before been observed on similar Hamiltonian 
models, in a noticeable way, for such a large range of energies. For example, 
for the XY (n = 2) version of the present model, such behavior was observed 
only near criticality. It is shown also that the (metastable state) maximum 
Lyapunov exponent decreases with N like A 

max ~ N K , where for the initial 
conditions employed herein (maximal magnetization), k = 0.225 ± 0.030, 
both above and below the critical point. 

Keywords: Hamiltonian dynamics; Heisenberg model; Long-range interac- 
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1. Introduction 



The Boltzmann-Gibbs (BG) statistical-mechanics represents one of the most successful 
physical theories, providing a good description of many experimental systems at equilib- 
rium [1-3]. The applicability of such a formalism is justified upon the validity of the 
ergodic hypothesis, which requires that the whole phase space should be equally visited 
in the infinite-time limit. Typically, this occurs for large Hamiltonian systems, with dy- 
namical variables connected by short-range interactions, leading to a microscopic dynam- 
ics characterized by a quasi-continuum Lyapunov spectrum, whose largest value remains 
positive in the thermodynamic limit. Such a chaotic dynamical behavior leads to a quick 
occupation of phase space, ensuring a safe use of the ensemble theory, and so, the BG 
statistics may be connected to the standard extensive thermodynamics in the well-known 
elegant manner. However, ergodicity and, consequently, the validity of the standard equi- 
librium ensembles, depend crucially on the nature of the Hamiltonian system considered 
[4]. In particular, systems characterized by long-range interactions, or long-range micro- 
scopic memory, may present a breakdown of ergodicity, leading to a fractal (or even more 
complex) structure in phase space. In such cases, the BG statistical-mechanics framework 
looses its validity and some more general theory must be employed. Recently, a large va- 
riety of evidences has been presented, exhibiting results that do not conform with the 
BG formalism; one may mention observations on turbulent plasmas [5], turbulent fluids 
[6-8], astrophysical systems [9-13], quantum chaos [14], logistic map [15], glasses [16,17], 
and complex systems [18,19], among others. Such systems exhibit evident inconsistencies 
with one of the main characteristics in the BG formalism, which is the extensivity of the 
entropy and the exponential weight factor associated with it. For an iV-particle system, 
the extensivity property means that quantities like the internal or free energy per par- 
ticle should approach a well-defined finite thermodynamic limit when N — > oo. When 
long-range interactions are present, each microscopic constituent of the system interacts 
with all the others, leading to an energy that depends more than linearly with N, and, 
obviously, to a nonextensive behavior. As a consequence of this, additivity does not hold, 
in the sense that the free energy of the whole system is not equal to the sum of the free 
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energies of its macroscopic parts, and the application of the BG formalism becomes ques- 
tionable. A generalized nonextensive thermostatistical formalism, proposed over a decade 
ago [20,21], seems to be a good candidate to deal with such systems. 

Recently, a lot of attention has been dedicated to a classical Hamiltonian system, 
namely, the inertial long-range-interaction XY model, which consists in an assembly of N 
classical planar rotators interacting through a long-range potential [22-35]. Such a system, 
which has been investigated numerically within the microcanonical ensemble, presented 
clear indications that a more general statistical-mechanics formalism is required for its 
description. In particular, for the case of infinite-range interactions, i.e., mean-field limit, 
for which a well-known continuous phase transition occurs, if one considers a total energy 
close to and below the critical energy, there exists a basin of attraction for the initial 
conditions for which the system gets captured in a metastable state, whose duration 
increases with N, before attaining the terminal thermal equilibrium. Therefore, if one 
considers the thermodynamic limit (N — > oo) before the long-time limit, the system will 
remain in the metastable state and will never reach the terminal equilibrium state, in 
such a way that the phase space will not be equally and completely covered. Moreover, in 
such a metastable state, the maximum Lyapunov exponent approaches zero, as N — > oo, 
contrary to what is expected in a standard BG equilibrium state. 

In the present work we perform molecular dynamical investigations, of the iso- 
lated inertial infinite-range-interaction Heisenberg ferromagnet, defined by N classical 
Heisenberg-like rotators. We show that the metastable state that occurs in the corre- 
sponding XY version of the model, near criticality, now appears in a noticeable way for a 
much wider extent, for an energy range that starts right below criticality and prolongs up 
to very high energies. In the next section, we review the equilibrium canonical-ensemble 
solution of the model. In section 3 we present the results of our numerical investigation, 
and finally, in section 4, we present our main conclusions. 
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2. The Equilibrium Canonical-Ensemble Solution 



Let us now work out the equilibrium canonical-ensemble solution of the model. For 
the sake of generality, throughout this section, we will deal with an inertial n-vector 
ferromagnet, composed by N classical rotators, each of them defined in an n-dimensional 
configurational space. The Hamiltonian of the system is given by 

i N n i N 

1=1 ^=1 1,3=1 
i N n i AT / n \ 

where the index \i (/i — 1, 2, • • • , n) denotes Cartesian components and L ifJl represents the 
yU-component of the angular momentum (or the rotational velocity, since we are assuming 
unit inertial moments) of rotator i. It is important to remind that the iV-dependence 
in the coupling constant above, is usually (but not necessarily, see [24]) introduced in 
order to yield a sensible thermodynamic limit, i.e, a finite free energy per particle when 
iV — > oo, within the equilibrium ensemble theory of standard statistical mechanics. 

The rotators are allowed to vary their directions continuously inside an n-dimensional 
sphere of unit radius, leading to the constraint 

fX = l (< = 1,2,...JV). (2.2) 

It should be mentioned that such a constraint reduces the number of degrees of freedom 
per particle to n — 1, in such a way that the total number of degrees of freedom of the 
system is given by N(n — 1). The model defined above recovers, as particular cases, the 
mean-field inertial XY (whose equilibrium canonical-ensemble solution was presented in 
Ref. [22]) and Heisenberg (whose dynamics will be discussed in the next section) models, 
for n = 2 and n = 3, respectively. Although our model is composed by one-dimensional 
inertial constituents (rotators), we shall sometimes refer to S = (Si, S2, ■ ■ ■ , S n ) as spin 
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variables, considering the close analogy of the above- defined model with the standard 
n- vector ferromagnet. 

One may now follow the standard procedure by rewriting the Hamiltonian as 



N 1 N n I 



N 



2 2 



i=i /x=i 



ZiV fi=i \i=i 



(2.3) 



in such a way that the partition function becomes 



Z = exp 
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N n \ / N n-1 
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(2.4) 



The constraints of Eq. (2.2) are taken into account in the equation above, through the 
delta functions, as well as in the products over the angular-momentum variables, which 
apply only over the effective number of degrees of freedom per particle, /i — 1,2 • • • ,n — 1. 
The squared J2i ma Y be linearized through the application of a Hubbard-Stratonovich-like 
transformation [2], which introduces a set of parameters {x^}; rescaling {x^} \fN{x IJ } 
one gets 



/ 2 xJV(-i)/ 2 



n 



'-Mr?/ 



N 



X 



nj/n [ rf ^p o v x^)] s nr 5* - 1 

M =1 \M =1 



(2.5) 



As usual, the site index may be discarded, and straightforward calculations lead to 



Z = 



AT(n-l)/2 



exp 



-f3N s 



n 
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, (2.6) 



where 
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1/2 



£ = 2(- 2 )/ 2 7rV (B - 2)/2 /(n-2)/ 2 (y); y = ( /? £ 4 



(2.7) 



with 7fc(z/) denoting modified Bessel functions of the first kind of order k. Considering N 
large, one may use steepest descents, 



/ 27r XAT(„-l)/2 / 

Z^ij expl — lexp 



-Nm&Xy (|^ -m£(y) 



(2.8) 



with the condition of maximum leading to the self-consistent equation, 



\ZdyJ I{n-2)/2 



(v)' 



(2.9) 



In the thermodynamic limit, one obtains the free-energy per particle, 



(2.10) 



as well as the internal energy per particle, 



^ — 1 1/-, -*2\ 



(2.11) 



where m represents the magnetization per particle, whose modulus is directly related to 
the parameter y, 



m = \m\ — — — 



13 I( n -2)/2(/3m)' 



(2.12) 



The critical temperature of the model may be obtained by considering m small and 
expanding the r.h.s. of Eq. (2.9) in power series (we work in units of ks — 1), 
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1 r(n/2) 1 

c " 2T[(n + 2)/2] ~n' (2 ' 13) 



which may be substituted into Eq. (2.11) to yield the critical energy density 



u c =i[l + (n-l)T c ] = l-^ (2.14) 



The above results recover those obtainted in Ref. [22], for n — 2, e.g., the self-consistent 
equation, 

leading to T c = 1/2 and w c = 3/4. For the case n — 3, one gets the well-known self- 
consistent equation, 

777, = cotanh(/3m) — — !— , (2.16) 

pm 

with T c = 1/3 and u c = 5/6. For n increasing from unit (Ising model) to infinity (spherical 
model), T c decreases form 1 to 0, and u c increases from 1/2 to I. 

In the next section we discuss the dynamics of the particular case n = 3 of the model 
defined above. 



3. Dynamics of the Mean-Field Inertial Heisenberg Model 
3.1 - Molecular Dynamics 

In this section we will present the results obtained by simulations of the constant-energy 
dynamics of the model defined by Eqs. (2.1) and (2.2) for the particular case n — 3, 
i.e., Heisenberg-like rotators. From now on, we will use the standard notation for the 
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Cartesian components of a Heisenberg model, i.e., /j, — x,y, z. The results to be discussed 
below were obtained by a direct integration of the equations of motion, 



Li-SiX - 



(< = 1,2,..-, 



N), 



(3.1a) 




(i = l,2,...,JV), 



(3.1b) 



which correspond to a set of 6iV equations to be solved numerically. 

For solving such a set of equations we have used a fourth-order Runge-Kutta-Merson 
integrator [36] with a time step of 0.05, leading, respectively, to the relative energy and 
spin-normalization conservations of 10~ 4 and 10~ 3 , or better. The total initial kinetic 
energy was divided into three equal parts, each of them to be assigned to a given set 
of Cartesian components of angular velocities {L^} (i = 1,2, •■-,N). We have always 
started the system with the so-called water-bag initial conditions [29,30,34] for each set 
of components of angular velocities, i.e., each set {L ifJ } was extracted from a symmetric 
uniform distribution and then, translated and rescaled to have zero total momentum. In 
what concerns the spin variables, we have started our simulations with all spins aligned 
along the z-axis (zero initial potential energy). Our measured quantities correspond to 
averages over N s distinct samples, i.e., different initial sets of {L^}. 

3.2 - Caloric Curve and Metastability 

In Fig. 1 we exhibit the caloric curve (full line) obtained by solving the equilibrium 
canonical-ensemble equations [Eqs. (2.11) and (2.12)]. We have chosen four particular 
values of the energy density to investigate, within our microcanonical-ensemble molecular 
dynamical approach, how (K)/N evolves in time, for different values of iV (it should 
be mentioned that, the quantity {K)/N, when evaluated at equilibrium, is expected 
to coincide with the temperature). Two of the chosen energies, u = 0.75 [Fig. 2(a)] 
and u = 0.96 [Fig. 2(b)], correspond, respectively, to values slightly below and above 
the critical internal energy (u c = 5/6). The energy u = 1.32 [Fig. 2(c)] is inside a 
range of energies where the kinetic temperature (K)/N presents a maximum discrepancy 
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between its values at intermediate and long times. Finally, the fourth chosen energy, 
u = 7.2 [Fig. 2(d)] corresponds to a value far above u c . In all the plots of Figs. 2(a)- 
(d) the system was started with the above-mentioned initial conditions and we have 
considered N s = 16 (N = 200), N s = 12 (N = 400), N s = 8 (N = 800), and N s = 4 
(N = 1600). We have observed that, after a short transient, the system rapidly reached 
a metastable, or quasi-stationary state (QSS), with a value of (K)/N higher than the 
one predicted by the canonical-ensemble equilibrium theory. It is important to remind 
that a QSS has been found also for the corresponding XY version of the present model, 
in an unambiguous way, only near criticality [23,25,27,29,30,34]. Except for very low 
energies, the QSS is always detected easily in the Heisenberg case; this is shown in Fig. 
1, where we exhibit the values of (K)/N for both QSSs (empty squares) and terminal 
equilibrium states (empty circles), for systems with N = 400 and different values of the 
energy density. One clearly sees that, in what concerns the value of (K) /N at the terminal 
equilibrium state, the thermodynamic limit may be attained within our computational 
effort (in the sense that the present microcanonical-ensemble numerical approach agrees 
with the equilibrium canonical-ensemble results) for u ^> u c , whereas for u ~ u c , one 
observes strong finite-size effects. In fact, the terminal-equilibrium value of (K)/N seems 
to reach its thermodynamic limit for very small systems, if u ^> u c [as shown in Fig. 2(d), 
for the case u = 7.2, where the size N = 200 appears to have converged], whereas near 
criticality (u ~ u c ) one may observe fluctuations [as shown in Figs. 2(a) and 2(b)]. 

Although large fluctuations may be observed on the values of (K)/N in the QSS, it 
appears evident that the gap with respect to the corresponding terminal-equilibrium-state 
values survives in the thermodynamic limit. If one defines the lifetime of the QSS (£qss) 
as the time at which (K)/N presents its half-way between the values at the QSS and the 
terminal thermal equilibrium, one concludes that such a quantity increases, essentially, 
linearly with N, as shown in Fig. 3, for the energies considered in Figs. 2(a)-(d). There- 
fore, the duration of the QSS increases with N, in such a way that, if the thermodynamic 
limit is performed before the long-time limit, the system will never relax to the terminal 
thermal equilibrium. 
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3.3 - Sensitivity to the Initial Conditions 

Let us now investigate how the maximal Lyapunov exponent, A max scales with N. Herein 
we shall use the well-known method for calculating such a quantity by considering the 
limit [37], 



where d(t) represents the metric distance calculated from infinitesimal displacements in 
phase space, at time t. We have carried simulations up to t = 10000 (i.e., 200000 time 
steps), for different system sizes (N = 50,100,400,1600,3200), and A max was obtained 
after averaging over N s = 50 samples for the smallest size (N = 50), whereas for all 
other sizes we have considered N s = 10. For the cases where there is an apparent QSS, 
the time interval considered ensures that the A max computed does indeed correspond to 
a quantity in the QSS [see, e.g., Figs. 2(a)-(d)], whereas for the cases where there is no 
evident QSS, the time used is expected to be sufficient for the system to have reached its 
terminal thermal equilibrium. As shown in Fig. 4, one has that A max ~ N~ K , similarly 
to what happens for the XY version of the present model [23,30,34,35]. We computed 
k for energies below u c (u = 0.30), as well as for u 3> u c (u = 7.2), with the above- 
mentioned initial conditions; we have obtained essentially the same estimate in both 
cases, k = 0.225 ± 0.030. It is important to remind that, in the mean-field inertial XY 
model, one has k = 1/3 for « > « c [23,30,34], whereas k — 1/9 for u < u c [35]; in 
the former case, there is no apparent QSS, in such a way that the estimate k = 1/3 is 
expected to apply to a terminal thermal equilibrium state, whereas in the latter, k was 
obtained at the peculiar QSS. Such estimates find no similar with those presented herein 
for the Heisenberg case, which apply to the QSS for u ^> u c , whereas for the low energy 
considered (u = 0.30) in the case u < u c we have found no clear evidence of a QSS 
(although the presence of a QSS, with a small gap, indiscernible due to fluctuations in 
{K)/N, with respect to the terminal equilibrium state, is not ruled out). However, the 
common feature observed for both XY and Heisenberg models, A max ~ N~ K (yielding a 



A 



max 




lim A(f), 



(3.2a) 




(3.2b) 
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zero maximal Lyapunov exponent in the thermodynamic limit), which seems to hold in 
the presence or not of an evident QSS, does certainly contradict the standard BG theory 
for an equilibrium state (which requires a finite A max in the thermodynamic limit). 

4. Conclusion 

We have analysed a system of N Heisenberg-like classical rotators with ferromagnetic 
infinite-range interactions. The dynamics of the model was studied within the micro- 
canonical ensemble, by directly solving the equations of motion. For a finite N, the time 
evolution of the kinetic temperature shows that there is a basin of attraction for the initial 
conditions for which the system gets caught in a metastable state, before reaching the ter- 
minal thermal equilibrium. We have shown that the duration of such a metastable state 
diverges as N — > oo, tipically linearly with N. Therefore, if the thermodynamic limit is 
considered before the long-time limit, the system will never relax to the terminal thermal 
equilibrium. We have also calculated the maximum Lyapunov exponent, above and below 
the critical point; in both cases, the scaling, A max ~ N~ K , was verified. For the particular 
initial conditions considered (maximal magnetization), the exponent k presents the same 
value for energies chosen above and below the critical point, k = 0.225 ±0.030. Above the 
critical point our estimate applies to the metastable state, whereas the estimate below 
the critical point is expected to hold for the terminal thermal equilibrium, since in such a 
case we have found no clear evidence of the existence of a metastable state. Preliminary 
studies suggest that, above criticality, the exponent k does not depend, within the error 
bars, on the initial conditions employed; however, below criticality, different initial condi- 
tions may possibly lead to a breakdown of universality, with different estimates for k. In 
particular, the k estimates below criticality seem to vary if the initial conditions for the 
spin variables break or not the Heisenberg-like symmetry of the system, i.e., if we start 
with m/0 (present paper) or m = 0. Further studies to clarify this point constitute the 
next step along the present lines. 
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Figure Captions 



Fig. 1: The caloric curve for the inertial mean- field ferromagnetic Heisenberg model 
obtained by the equilibrium canonical-ensemble solution (full line). In the vertical axis, 
(K)/N, is a quantity that is expected to coincide with the temperature, at equilibrium. 
The dashed vertical line signals the second-order phase transition critical energy density, 
u c = 5/6. The empty squares and circles represent the estimates of (K)/N from the 
microcanonical-ensemble numerical analysis of a system of size N = 400, at the metastable 
state and after that, respectively. 

Fig. 2: The microcanonical time evolution of (K)/N is represented for several 
system sizes and different energy densities: just below criticality [u = 0.75 (a)], just 
above criticality [u = 0.96 (b)], at the region where the gap between the metastable and 
terminal equilibrium states is maximum [u = 1.32 (c)] and for u >> u c [u — 7.2 (d)]. The 
initial conditions are water-bag for velocities and m—1. 

Fig. 3: Log-log plots of the lifetime (£qss) of the QSS as a function of N, for the 
energy densities considered in Figs. 2(a)-(d). Simple linear fits yield the slopes 0.99±0.05 
( u = 0.75), 1.10 ± 0.06 (u = 0.96), 1.01 ± 0.02 (u = 1.32), and 0.95 ± 0.06 (u = 7.20). In 
each case, the slope is very close to 1 (represented by the dashed line), in such a way that 
£qss ~ N in all cases. The initial conditions are water-bag for velocities and m — 1. 

Fig. 4: Log-log plots of the maximum Lyapunov exponent A max versus N, for 

energies densities below (u = 0.3) and above the critical point (u = 7.2), showing the 
scaling A max ~ N~ K . The slopes are essentially the same, within the error bars, yielding 
k = 0.225 ± 0.030. For u = 7.2, the distinction between the metastable and teminal 
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equilibrium states is very clear, and this value of k corresponds to the metastable one. 
For u = 0.3, the distinction in not very clear, and this value of k presumably corresponds 
to the terminal equilibrium state. 
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